A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
02 (12.03.4 12.12.1 03.13.6) 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



On 
On 



o 



o 
o 



Of 
I 1 

o 



>< 



The accuracy of parameters determined with the 
core-sampling method: application to Voronoi tessellations 

Andrei G. Doroshkevich 1 , Stefan Gottlober 2 , and S0ren Madsen 3 

Theoretical Astrophysics Center, Juliane Maries Vej 30, DK-2100 Copenhagen 0, Denmark 
Astrophysikalisches Institut Potsdam, An der Sternwarte 16, D- 14482 Potsdam, Germany 
3 Copenhagen University, Astronomical Observatory, Juliane Maries Vej 30, DK-2100 Copenhagen 0, Denmark 

Received ; accepted 



Abstract. The large-scale matter distribution represents 
a complex network of structure elements such as voids, 
clusters, filaments, and sheets. This network is spanned 
by a point distribution. The global properties of the point 
process can be measured by different statistical methods, 
which, however, do not describe directly the structure el- 
ements. 

The morphology of structure elements is an impor- 
tant property of the point distribution. Here we apply 
the core-sampling method to various Voronoi tessella- 
tions. Using the core-sampling method we identify one- 
and two-dimensional structure elements (filaments and 
sheets) in these Voronoi tessellations and reconstruct their 
mean separation along random straight lines. We com- 
pare the results of the core-sampling method with the a 
priori known structure elements of the Voronoi tessella- 
tions under consideration and find good agreement be- 
tween the expected and found structure parameters, even 
in the presence of substantial noise. We conclude that the 
core-sampling method is a potentially powerful tool to in- 
vestigate the distribution of such structure elements like 
filaments and walls of galaxies. 



1. Introduction 

The first deep galaxy surveys which have become avail- 
able have shown that the galaxies are not at all homo- 
geneously distributed on large scales. Large voids, small 
filaments, massive clusters, and walls of galaxies can be 
observed. How these structures have evolved from the ini- 
tially nearly homogeneous universe is one of the central 
problems in modern cosmology. Many theoretical mod- 
els of structure formation have been suggested. Powerful 
mathematical statistical methods are necessary in order 
to describe this structure, to obtain its quantitative char- 
acteristics, to compare the theoretical models and simu- 



lations with observational catalogues, and, eventually, to 
find the correct theory of structure formation. 

The standard approach for testing models is to define 
a point process which can be characterized by statisti- 
cal methods. This could be the distribution of galaxies 
of a specific type in deep surveys or clusters of galaxies. 
In order to compare models of structure formation, the 
different distribution of dark matter particles in N-body 
simulations could be tested as well. 

The most widely used statistics are the n-point cor- 
relation functions, the counts-in-cells method, and the 
void probability function. The geometrical and topological 
properties of the point distribution can be investigated by 
the percolation technique (Zeldovich, Einasto, Shandarin 
1982), the minimal spanning tree (Barrow, Bhavsar, & 
Sonoda 1985), the genus of the smoothed density field 
(Gott et al. 1986, 1989), and the Minkowski functionals 
(Mecke et al. 1994). 

All these different statistics measure the global prop- 
erties of a point process. However, when looking at a point 
process we also see different structures. The morphology 
of these structures is an important property of the point 
distribution. (Note also the warning of Barrow & Bhavsar 
(1987) that we tend to see by eye structures which do 
not exist). Therefore, it is important to develop statistics 
which distinguish between different structure elements in 
a point distribution and find their typical scales (if any 
exists). 

In this direction a first attempt was made by Vish- 
niac (1986). He used the moments of a point distribution 
in a window to measure the amount of filaments in two- 
dimensional galaxy distributions. Recently, this method 
was generalized to three-dimensional point distributions 
(Luo, Vishniac 1995). Both the number of filaments and 
walls and typical scales of the structure can be found. 
Minkowski functionals are also efficient discriminators for 
idealized one-, two- or three-dimensional structure ele- 
ments (Schmalzing, Kerscher, Buchert 1996). However, 
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using Minkowski junctionals it is difficult to extract such 
structure elements from a superposition of all possible el- 
ements including Poisson noise. 

Using the core-sampling method (Buryak, Demiahski, 
& Doroshkevich 1991), Buryak, Doroshkevich & Fong 
1994) one can find structure elements and their typi- 
cal separation. The method was designed to find one- 
dimensional filaments and two-dimensional walls within 
observational surveys of galaxies (Buryak, Doroshkevich 
& Fong, 1994, Doroshkevich et al. 1996b) and in simulated 
samples (Doroshkevich et al. 1996a). Here we want to ap- 
ply the method to an idealized model mimicking many 
features of the observed patterns in the galaxy distribu- 
tion. The model employed here is based on the concept 
of a Voronoi tessellation (Goldwirth, da Costa, & Van de 
Weygaert (1995)). Such a test allows us to compare the 
input and output structure parameters and, thus, to test 
how powerfully this method can discriminate structure el- 
ements, determine their distribution and typical parame- 
ters. 

The basic idea of the core-sampling method is to re- 
duce the analysis of a three-dimensional galaxy distribu- 
tion to the investigation of the distribution of structure 
elements along random straight lines (the 'core'). Such 
an approach allows us to avoid any discussions and de- 
scriptions of the very complicated multiconnected struc- 
ture as a whole as well as any conventional definition of 
a 'structure element' and 'void' in the three-dimensional 
space. The core-sampling method has to deal only with a 
one-dimensional point distribution. It allows us to define 
two typical populations of structure elements, namely, fila- 
ments and sheets, and provides us with two objective char- 
acteristics of the spatial distribution of these elements, 
namely the surface density of filaments and the linear den- 
sity of sheets. The surface density of filaments Uf is simply 
the mean number of filaments intersecting a unit area of 
arbitrary orientation, while the linear density of sheets 
o~ s is simply the mean number density of sheets crossing 
an arbitrary straight line. These characteristics provide a 
local description of a random three-dimensional network 
structure. In Sect. 2 we will briefly outline the general 
properties of this method 

The three-dimensional Voronoi tessellation consists of 
a number of Voronoi cells enclosed by the planes bisect- 
ing the lines between the nuclei of neighbouring cells (van 
de Weygaert, 1991). The intersections of these planes are 
lines and points. Therefore, the three-dimensional Voronoi 
tessellations consists of three structure elements: the walls 
around the voids, the edges (intersection of walls) and the 
nodes. Within the three-dimensional Voronoi tessellation 
particles may be distributed on each of these structure el- 
ements or distributed around them according to a given 
distribution function. Moreover, the Voronoi tessellation 
can be superimposed by randomly distributed particles. 
Thus it is a unique tool for testing algorithms purporting 
to describe the structure elements of a point distribution. 



In Sect. 3 we will briefly describe Voronoi tessellations 
and the creation of synthetic samples. We apply the core- 
sampling method to these synthetic samples and compare 
the results with the input structure element. We discuss 
the accordance between the expected and found structure 
parameters and conclude that the core-sampling method is 
a powerful tool to investigate the distribution of structure 
elements. 

2. The Core-Sampling Method 

2.1. Basis of core- sampling method 

The mathematical basis for this method had been de- 
scribed by Buryak, Demiahski, & Doroshkevich (1991) and 
Buryak, Doroshkevich & Fong (1994). A recent detailed 
discussion can be found in Doroshkevich et al. (1996b). 

First let us give the salient points of the core-sampling 
method. 

1. The distribution of structure elements along a random 
straight line is assumed to be Poissonian. Thus, a ID 
cluster analysis is utilized to discriminate the structure 
elements among the sample of points and to find the 
number and the mean separation of structure elements. 

2. For the following cluster analysis the fields are then all 
organized into an 'equivalent single field' by combining 
the separate ID distributions one after the other along 
a line, with the first point of a field placed on top of 
the last point of the preceding field. 

3. The dependence of the number and mean separation 
of structure elements on the diameter of the core al- 
lows a rough discrimination between the filament and 
sheet-like populations of the structure elements and 
yields the fundamental characteristics of the structure, 
namely, the surface density of filaments, 07, and the 
linear density of sheets, er s . 

4. The sample of points under consideration is reduced 
by rejecting poorer and sparser structure elements. In 
this manner, the mean characteristics can be found as 
a function of the threshold richness. 

Let us emphasize here that in practice the one- 
dimcnsional analysis is very convenient in many respects. 
We are using here and in numerical simulations cylinders 
around straight lines, but they can be easily replaced by 
cones for observational surveys. Thus, the core-sampling 
method can be directly used to analyze pencil beam sur- 
veys as well as almost two-dimensional samples like the 
slices of the deep Las Campanas survey or real three- 
dimensional surveys. Moreover, in case of two- and three- 
dimensional samples it allows measurement of structure 
parameters in different directions, for instance, both along 
the line of sight and the transverse direction. Thus, red- 
shift space distortions can be extracted (Doroshkevich et 
al. 1996b). 

The assumption of a functional form for the distri- 
bution of structure elements along the core is essential for 
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our method, because it enables transformation of the point 
distribution into a distribution of structure elements. Once 
this transformation is accomplished, the appropriateness 
of the assumed functional form can be tested and the val- 
ues of the functional parameters can be determined. 

As stated above, we assume a ID Poissonian law for 
the distribution of structure elements (not galaxies) along 
the axis of cylinder. It is a simple distribution, and we 
shall see that it is also a reasonable assumption. Indeed, 
the Poissonian distribution arises naturally for some theo- 
retical models (White 1979, Buryak et al. 1991) when the 
mean separation of structure element exceeds the corre- 
lation length. The validity of this assumption, however, 
cannot be tested a priori, and a possible difference be- 
tween the assumed and the actual distributions depends 
on the sample in question, thus limiting the precision of 
the final results. This assumption was valid for all cases 
in previous investigations. Here we are testing the method 
with Voronoi tessellations, for which our assumption has 
to be tested as well. 

2.2. Sample preparation 

The first step of core-sampling analysis is the prepara- 
tion of a sample with suitable parameters. This sample 
depends both on the radius of the cylinder R c and the 
linking length l c of the culling procedure. The sample is 
the basis for the further analysis. 

The diameter of the cylinder must be a few times 
smaller than the expected size of the cells in order to avoid 
a masking effect due to possible overlapping of projected 
elements. On the other hand, the cylinder must be wide 
enough to allow during further analysis the sequential re- 
duction of the radius by a factor of at least 2, prefer- 
ably more, because the radius is that diagnostic parame- 
ter which allows one to distinguish between filaments and 
sheets. The number of particles has to be large enough to 
guarantee stability and reliability of results even for the 
smallest radius used for analysis. In practice one observes 
that this can be achieved if a structure element contains, 
in the mean, at least 5, or preferably more, particles. 

The final step of the sample preparation is sample re- 
duction, which is performed sequentially during the anal- 
ysis. Sample reduction provides an additional method of 
distinguishing between different populations of filaments 
and sheets according to the mean density or, in other 
words, according to the multiplicity of the structure el- 
ements. It has to be performed before any further steps of 
analysis. 

To do the sample reduction, a ID cluster analysis of 
each field is performed with decreasing linking length l c 
of the culling procedure. All clusters with fewer mem- 
bers than a constant threshold multiplicity fx t hr are re- 
jected. Then the remaining number of particles in the sam- 
ple N r depends on the current linking length l c and the 
threshold [i t hr ■ For each linking length l c a certain fraction 
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/ = N r /N a of points remains (N a is the total number of 
particles). It is convenient to use the fraction / together 
with the threshold multiplicity fithr, as parameters char- 
acterizing the reduced sample. 

This procedure reduces the full sample, rejecting more 
and more points associated with poorer or sparse clumps. 
Only the tighter clumps are retained for further analy- 
sis. This approach emphasizes the local density within 
the structure elements. It allows the rejection of the low- 
density haloes from the structure elements and sparse 
structure elements as a whole. 

2.3. Detection of the mean separation of structure ele- 
ments 

The next step of analysis is the detection of the mean 
separation of structure elements. At first, the axes of all 
cylinders are randomly combined to one line by identify- 
ing the last point of one axis with the first of the next. 
At this line a ID cluster analysis with increasing linking 
length R is performed. Assuming a Poisson distribution, 
the number of clusters Nr as a function of R is given by 
the relation 

ln(N R ) = ln(A ) - R/R . (1) 

For a true Poissonian sample, these values are related to 
the length of the sample, Do, defined by the nearest and 
the farthest points in an average field, via the relationship 



No/D = 1/R . (2) 

Thus, the difference between the linear density of struc- 
ture elements Nq/Dq and their mean separation Rq ob- 
tained from Eqs. (|l|) and (||) can be considered a measure 
of the systematic error due to the deviation of the actual 
from the assumed (Poissonian) distribution of structure 
elements along the analyzed cylinder. To decrease this er- 
ror we use an automatic procedure to find the optimal 
interval for R for the fit to Eq. (|l|). The upper limit to 
this interval is fixed by the input of a minimum number 
Nr > N m i n , whereas the lower limit can vary and is de- 
fined by the condition 

1 - e < R N /D <l + e (3) 

where e is a desirable precision of fitting parameters. In 
general, a reasonable precision (~10%) can be achieved in 
the analysis. 

If a desirable precision cannot be achieved, the reason 
is usually that the distribution used is far from Poissonian 
and the sample under consideration must be changed (e.g., 
by changing the cylinder radius). 
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2.4- Identification of structure elements 

The final step of analysis is the discrimination of filaments 
and sheet-like structure elements and determination of the 
parameters cr/ and a s for both populations. To this end, 
we use a simple geometrical model for the structure ele- 
ments (Buryak, Doroshkevich & Fong (1994)). According 
to this model, in each narrow cylinder the structure can 
be considered as a system of randomly distributed lines 
(filaments) and planes (sheets) which contain all points. 
In this model, the filaments are considered straight lines, 
and the sheets are considered flat planes. Of course, this 
approach is limited, and it cannot be used as an accu- 
rate description of the true matter distribution on large 
scales; one can best characterize it as an intermediate step 
between the local description of the matter distribution 
with density and velocity fields and the global description 
obtained by the topology or Minkowski techniques. 

As it was mentioned above, we characterize the ran- 
dom distribution of straight lines (filaments) by their sur- 
face density, cr/, i.e., the mean number of lines intersecting 
an unit area of arbitrary orientation. The random distri- 
bution of planes can be characterized by the linear density, 
a s , i.e., the mean number density of planes (sheets) cross- 
ing an arbitrary straight line. Equivalently, we can use the 
values D s = a J 1 and Df — a^ as typical measures of 
the mean separation of structure elements. To character- 
ize the structure as a whole we can estimate the mean 
separation of structure elements — filaments and sheets 
combined — by < Df s >, which can be thought of as the 
diameter of a sphere containing, on average, two structure 
elements, 



The main characteristics of the structure can be ob- 
tained by fits to the radial dependence of the linear den- 
sity of clusters, No/ Do, and to the radial dependence of 
the mean separation of structure elements, i?o, 

N /D =a s +<n Rcaf, (5) 
and 

1/Ro = er s + 7T R c o-f, (6) 

over the range of variation of core radius R c . Here Dq is 
the mean depth of field. 

Clearly, these parameters also depend on the sample 
being studied. In our analysis, Eqs. @ and (||) were fit- 
ted by a maximum likelihood technique, and the resulting 
mean values of at and o~ s were accepted as the final esti- 
mates of the structure parameters. The difference of the 
estimates from Eqs. (|j) and (Q) was included into the er- 
rors of the final values. 



2.5. Methodological remarks 

For the idealized model considered above, the core- 
sampling method would be expected to reproduce well the 
characteristics of the sample under consideration. How- 
ever, in reality, several factors may distort the final re- 
sults. For example, filaments may be apparently detected 
as sheets, if they pass through the center of the core. Also 
the noise of randomly distributed particles must be taken 
into account. Therefore, a specific technique must be used 
to obtain stable results. 

The analysis of the Las Campanas Redshift Survey 
(Doroshkevich et al. 1996b) has shown that the stability 
of core-sampling with respect to various accidental or sys- 
tematic variations of density is very high. During the se- 
quential random rejection of particles, the final estimates 
of mean separation of structure elements increased only 
as / -1 / 3 . In general, the precision and stability of the 
core-sampling method depends on the density contrast in 
structure elements relative to the mean number density 
of points. In case of pure Poissonian point distribution 
the 'structure' parameters have been found with errors of 
about 50% (Buryak, Doroshkevich & Fong, 1994). Similar 
errors were found for the Soneira-Peebles model. On the 
contrary for the observed Las Campanas Redshift Survey 
errors were found to be about 10% only. 

1. Sample preparation. 

During the first step of analysis, when the sample is 
prepared, one has to check whether the resulting ID 
cluster distribution is Poissonian as assumed. One has 
to ensure this requirement by the variation of the ra- 
dius of cylinder and the range of fit in Eq. (|l|). During 
the further analysis the range of core radii used for the 
fit of Eqs. (||) and (||) is important to get the two val- 
ues in agreement. These questions have to be solved 
before the final analysis. This means that during each 
of the intermediate steps leading to Figs, [j] to |], the va- 
lidity of all assumptions must be checked. Otherwise, 
the final result will be meaningless. 

2. The surface density of filaments. 

The richness of a filament in the core depends both on 
the actual properties of the point distribution analyzed 
and the geometry, i.e. the position and orientation of 
the filament relative to the core. Only the first piece of 
information is of interest; the geometry can be charac- 
terized analytically and used to improve the stability 
of the method. 

The surface density of filaments af{f,(J*thr) depends 
on the fraction of particles retained in the sample and 
the multiplicity threshold. The main characteristic of 
filamentary structure is the full surface density of fil- 
aments, cr/(l,l). However, at this point the noise is 
usually high and it is desirable to improve the esti- 
mate of this value. 

Basing on the geometrical model described above one 
can calculate the distribution function for the length of 
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intersection of identical filaments with the core. Thus, 
calculating the dependence of both the fraction / and 
the surface density 07 on this length, we find in linear 
approximation 

^/(/) « - yffl - /) 1/2 + -■)■ (7) 

This relation estimates the fraction of poor filaments 
for the ideal case without noise. The linear fit of Eq. 

over some range of / allows one also to improve the 
estimate of 07 (1) for a moderate noisy sample. How- 
ever, in the case of very noise samples with a signifi- 
cant fraction of filaments consisting only of one point, 
it can be only used for small linking lengths, l c , of the 
culling procedure. Indeed, if l c < Rq the main part 
of noise particles will be rejected (see Fig. 8, right). 
Evidently, for smaller values of /, the function 
resembles the actual properties of the filaments (see, 
e.g. Doroshkevich et al. 1996b). 
3. The linear density of sheets. 

The richness of the elements of the sheet-like popu- 
lation both in observational data and simulations is 
usually quite high. Thus, the linear density of sheets 
can be estimated for a wide range of parameters / and 
Hthn providing more reliable results. However, even 
in this case the final results are sometimes distorted 
because during reduction the fields are shortened as a 
result of the appearance of empty edges. To avoid these 
distortions, a correction procedure (see Doroshkevich 
et al. (1996b) was used. 

The next problem is a filament population masquerad- 
ing as sheets. These apparent sheets are detected by 
the program if filaments cross the cylinder near its 
axes. Obviously, this effect will be proportional to the 
thickness of the filaments. The surface density cr" of 
these apparent sheets is a a s — nafrj, where 77 is the 
radius of the filament (Buryak, Doroshkevich & Fong, 
1994). Usually the influence of this effect is negligible, 
but it can become important for thick filaments. More- 
over, apparent sheets can be generated by the branch 
points (knots) of filamentary structure or other clumps 
in the point distribution if their size is comparable to 
the diameter of the core. The contribution of these 
knots can be estimated as cr" = Trn c i m rl lm , where n c i m 
and r c i m are respectively the volume density and the 
effective radius of the clump. For a random network 

structure, the volume density of the knots is approxi- 

3/2 

mately n chn » oy . 

3. Analysis of Voronoi- Tessellations 

3.1. Construction of Voronoi tessellations 

We have used the code developed by R. Van de Weygaert 
(1991, 1994, see also Goldwirth, da Costa, & Van de Wey- 
gaert (1995)) in order to construct Voronoi tessellations. 
Although the code has been developed to allow for much 



more sophisticated cases, we limit ourselves to the sim- 
plest case where a given number of points is distributed 
in walls or filaments of a given thickness. To mimic noise 
we have added an Poisson distribution of particles in the 
whole box. 

Using the same nuclei of the Voronoi tessellation, the 
filaments (defined as the edges of the walls) would have ap- 
proximately the same separation as the walls (contrary to 
what we expect in nature). Therefore, we did not combine 
a filamentary structure with a wall-like structure in one 
realization, but made a second realization for filaments 
alone. Note also that the filaments of a Voronoi tessel- 
lation are not distributed like filaments of real galaxies. 
Nevertheless the points of these Voronoi tessellation form 
one-dimensional structure elements with a priori known 
properties which should be recovered by the core sampling 
method. 

The mean diameter of the Voronoi cells and the size 
of the box are free input parameters of the code. For the 
core-sampling analysis it is only important that the mean 
number of structure elements along a random line crossing 
the box be large enough. We have chosen 6 cells and 20 
cells for the analysis of walls and filaments, respectively. 
The analysis itself can be performed in units of the box 
length. For the graphical representation we have chosen 
formally a box size of 300 Mpc. We have chosen 200000 
and 400000 points to represent the walls and the filaments, 
respectively. Due to the formal definition of the box size 
these numbers are not directly comparable with the num- 
ber of galaxies in the corresponding volume. However, in 
that case the mean number density of points is of the order 
10 -2 Mpc -3 , i.e of the same order as the galaxy density. 
Also the the mean separation of structure elements is of 
the same order as found for galaxies, so that we can in- 
deed draw conclusions from our test to the situation found 
in observational surveys. The noise level chosen (50% for 
walls and 25% for filaments) is in these realisations an ap- 
proximate upper limit for which the known properties can 
be recovered with high accuracy. 

3.2. Wall distribution 

For the first test, we chose a mean size of the cells of 
50 Mpc. Our first realization consists of 6 3 randomly dis- 
tributed cells in box of 300 Mpc length. 200 000 points 
are distributed only in walls around the cells. We chose a 
Gaussian density profile of the walls with a thickness of 
2 Mpc. On the left hand side of Fig. [jj a slice of 8 Mpc 
thickness of this realization is shown. In order to test the 
stability of the core-sampling method, we made a second 
realization with 100 000 additional points which were ran- 
domly distributed (Fig. [l], right). 

As explained in Sect. 2.1, the starting point of the 
core-sampling analysis is the construction of randomly dis- 
tributed cylinders in the simulation. For illustration, we 
put four of these cylinders into each of the slices shown 
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Fig. 1. Top: a slice of 8 Mpc thickness, bottom: radial particle distribution in the four cylinders shown at the top part of the 
figure (the dashed line indicates the minimum core radius used in the analysis), left: only walls, right: walls with a background 
of randomly distributed particles 



in Fig. [y. The diameter of the cylinders is equal to the 
thickness of the slice (8 Mpc). In the lower part of the fig- 
ures, the distribution of points in these cylinders is plotted 
vs. the radius of the cylinder. This radius will be used in 
the following as one diagnostic parameter as described in 
Sect. 2. Indeed, for sheet-like structures like the walls in 
our Voronoi tessellations, the mean number of clusters in 
the ID cluster analysis is independent of the core radius. 
One can clearly see from Fig. |l| that in both cases almost 
all walls will be identified as walls. 

Next, all particles N a in a given cylinder of radius R c 
are projected onto the axis of the cylinder. From Fig. [I] 
it is clear that there is also noise besides real structures 
in the resulting one-dimensional point distribution. The 
degree of noise depends on the chosen parameters of the 
generated sheets and on the parameters of the core. The 
noise particles in underdense regions can be removed us- 
ing the reduction procedure described in Sect. 2.1. The 
reduced sample at which the further analysis will be per- 
formed is characterized by three parameters: the fraction 
/, the threshold fithr, and the core radius R c . 

Further, we use the iterative fitting procedure de- 
scribed in Sect. 2.2 to find the best fit to Eq. (Q). Ac- 
cording to our assumption, the resulting clusters must be 
Poisson distributed for the some range of linking length. 
However, this assumption must be tested for the samples 
under consideration. As an example we show in Fig. || the 
fit to Eq. (|l]) for four sets of parameters (/ = 1, 0.5, r = 4 
Mpc, 2 Mpc and fi t hr = 2: left for walls and right for walls 
with randomly distributed points). 

Using Eqs. (||) and (m, where R c is the diagnostic pa- 
rameter, we determine now the surface density of filaments 
and the linear density of walls for each pair of parameters 
/ and Hthr- This procedure is illustrated by Fig. ||. for four 
pairs of parameters. (Note that each curve of Fig. || corre- 
sponds to one point in Fig. |^.) From Fig. ^| we extract the 
linear density of walls as the zero point of the curve and 
the surface density of filaments as the gradient. The fits 
to Eqs. (§) and (§ were performed over 9 equally spaced 
values of R c in the interval 2 Mpc < R c < 4 Mpc. 

As the final result of the analysis, we present in Fig. |3] 
the mean distance of the walls D s = aj 1 and in Fig. || the 
surface density of filaments aj depending on the threshold 
IHhr- (Note, that each curve of Fig. || provides one point 
in Fig. [| and one point in Fig. ||.) The final values of a s 
and it/ are listed in Table |l|. 

The noise of apparent filaments is an objective char- 
acteristic of the sample under consideration. The popu- 



lation of filaments dominates for all thresholds [ithr that 
resemble the small density contrast in the walls, since the 
number of particles concentrated into sheets exceeds the 
noise particles only by a factor of two. In this respect our 
synthetic models differ from the Las Campanas Redshift 
Survey, where for high threshold multiplicity the filament 
population had disappeared, and only the richest sheet- 
like elements survived. 

3.3. Filament distribution 

As a second test we have investigated a pure filamen- 
tary structure and filaments superimposed by randomly 
distributed particles. To this end we have constructed a 
Voronoi tessellation with 400 000 particles which are dis- 
tributed at the edges of the cells, and a second one with 
an additional 100 000 randomly distributed particles. For 
this test we chose a mean size of the cells of 15 Mpc. In 
Fig. H we show a slice of 4 Mpc thickness and four cores of 
radius 2 Mpc. Further steps in the analysis are as above. 
We made two realisations of Voronoi tessellations contain- 
ing only filamentary structure: In the first realisation we 
chose a radius of filaments of 0.2 Mpc, in the second 0.1 
Mpc. 

For these synthetic models, we omit the figures corre- 
sponding to Figs. HH but present only the final results, i.e. 
the density of walls cr s (f) and filaments o-f(f) vs. (1 — f) 1 ^ 2 
for two thresholds [i t hr = 2, 4 and the two different radii 
of filaments, (see Figs. j7f||). The final structure parame- 
ters are given in Table |f 

Let us note that in this case the mean edge separation 
along a random straight line is not identical to the mean 
size of the cell but is of the order of the mean length of 
the filaments. The mean distance of particles along the fil- 
aments in our sample is of the order of 1 Mpc. Therefore, 
we expect that the set of particles forms a broken net- 
work structure which can be characterized by the surface 
density of filaments cr/(/)- With a higher mean density 
we could create also a complete network structure which, 
however, in reality does not exist. 

Since all filaments were created using the same pro- 
cedure their parameters can vary only statistically. We 
can use Eq. ^ to determine 0/(1) as is shown in Fig. || 
(dashed line). In case of /ithr = 4 the points corresponding 
to a large fraction / do not follow the line predicted by 
Eq. (0). Rejecting these points, a linear fit would predict 
a somewhat smaller 07. We expect this result, because we 
are now measuring only the rare high-density filaments. 
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The noise strongly increases the full surface density 
of filaments <r/(l) (Fig. ||, right). The rejection of (10 - 
15)% of particles depresses the noise impact and leads us 
back to noiseless data. Since most or all observed galaxies 
are incorporated into structure elements the noise impact 
is not expected to be severe in catalogues. Following the 
rejection procedure we can obtain for all realizations the 
same result 0/(1) = 0.0054 which corresponds to 13.6 Mpc 
mean separation (see Fig. ||) . 

Fig. H (bottom) indicates that some of the filaments 
masquerade as sheets. We have discussed this problem al- 
ready in Sect. 2.4. In this case the error bars strongly 
increase as a result of the sparse sample statistics of el- 
ements which the core-sampling method wrongly detects 
as sheets (see Fig. fjj). 

4. Conclusions 

We have constructed a series of Voronoi tessellations with 
different input parameters. Our realisations consist either 
only of structure elements or of structure elements and 
noise. We have chosen a number of points comparable with 
the situation in observational samples and simulations of 
large scale structure formation, i.e. our realisations cor- 
respond to models with all galaxies within structure ele- 
ments and models with additional randomly distributed 
galaxies. 

In the first two tessellations, only walls were con- 
structed. Using the core-sampling method, we detected 
walls with a mean separation of 50 Mpc (see Table |l| and 
Fig. |2|) in these synthetic models. This corresponds exactly 
to the input parameter (the mean void size) of the Voronoi 
tesselation. The method provided an accurate reconstruc- 
tion despite the additional noise which was added to the 
tessellation (Fig.[| , right). On the other hand, the scales 
of the spurious filaments in the tessellation (Fig. ||) are 
sensitive to the addition of noise. The specific dependence 
found appears to be a property of the particular synthetic 
model generated. 

The next four tessellations are all characterized by a 
mean cell size of 15 Mpc. The cells are surrounded only by 
filaments. Two tessellations were made with filaments of 
radius 0.2 Mpc, two more with 0.1 Mpc. Independent of 
the radii of filaments, we were able to determine the same 
mean distance of filaments (13.6 Mpc); however, this esti- 
mate was slightly reduced by adding randomly distributed 
points as background (cp. Fig. || , dashed lines). The de- 
tected spurious walls have distances of the order of the 
box size and very large errors so that they can be easily 
ruled out. 

Our results demonstrate clearly that the core-sampling 
method allows one not only to determine the structure el- 
ements, but also to measure their characteristic mean sep- 
aration or density almost independently of the influence of 
a substantial noise component. From the good agreement 
between the expected and recovered structure parameters, 
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we conclude that the core-sampling method is a powerful 
tool for further investigations of observational surveys and 
the determination of structure elements such as filaments 
and walls of galaxies. Recently, the core-sampling method 
has been successfully applied both to numerical simula- 
tions (Doroshkevich et al. 1996a) and the Las Campanas 
galaxy sample (Doroshkevich et al. 1996b). 
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Table 1. Test of the wall distribution. 



thr &s O f a 3 Of 

lQ-^pc" 1 10~ 2 Mpc~ 2 lO^Mpc' 1 10~ 2 Mpc~ 2 

2 2.1 ±0.29 0.43 ±0.21 1.8 ± 0.20 2.7 ±0.16 

4 2.0 ± 0.13 0.43 ±0.19 2.5 ± 0.29 2.7 ±0.16 

6 1.8 ± 0.16 0.48 ±0.21 2.1 ±0.17 2.7 ±0.16 



Table 2. Test of the filament distribution. 



H-thr 


Us 

lO^Mpc- 1 


Of 

10^ 2 Mpc- 2 


noise 

lO^Mpc- 1 


noise 

f 

10^ 2 Mpc- 2 


r=0.2 Mpc 


2 


4.4 ±0.8 


0.57 ±0.02 


5.3 ± 1.0 


0.92 ± 0.05 


4 


7.3 ±0.8 


0.30 ±0.02 


6.5 ±0.6 


0.50 ± 0.04 


r=0.1 Mpc 


2 


4.5 ±0.9 


0.58 ±0.01 


3.5 ± 1.2 


0.90 ± 0.05 


4 


6.1 ±0.8 


0.31 ±0.01 


5.5 ± 1.0 


0.49 ± 0.04 
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Fig. 3. A fit to the Eqs. (5) and (6) from which No/Do) (dots) and 1/Ro (triangles) are obtained, (left: only walls, right: with 
a background) 
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Fig. 4. Mean separation of walls for three multiplicity thresholds (left: only walls, right: with background particles) 
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Fig. 5. Mean surface density of apparent filaments for three multiplicity thresholds (left:only walls, right: with background 
particles) 



Fig. 6. Top: a slice of 4 Mpc thickness, bottom: radial particle distribution in the four cylinders shown at the top part of the 
figure (the dashed line indicates the minimum core radius used in the analysis), left: only filaments, right: filaments with a 
background of randomly distributed particles 
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Fig. 7. Mean separation of apparent walls for two multiplicity thresholds (top: radius of filaments 0.2 Mpc, bottom: radius of 
filaments 0.1 Mpc, left: only filaments, right: with background particles) 
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Fig. 8. Mean surface density of filaments for two multiplicity thresholds (top: radius of filaments 0.2 Mpc, bottom: radius of 
filaments 0.1 Mpc, left: only filaments, right: with background particles) 



